Drying of complex suspensions 
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We investigate the 3D structure and drying dynamics of complex mixtures of emulsion droplets and colloidal 
particles, using confocal microscopy. Air invades and rapidly collapses large emulsion droplets, forcing their 
contents into the surrounding porous particle pack at a rate proportional to the square of the droplet radius. By 
contrast, small droplets do not collapse, but remain intact and are merely deformed. A simple model coupling 
the Laplace pressure to Darcy's law correctly estimates both the threshold radius separating these two behaviors, 
and the rate of large-droplet evacuation. Finally, we use these systems to make novel hierarchical structures. 
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The drying of suspensions of colloidal particles gives rise 
a plethora of fascinating phenomena, from the "coffee-ring" 
effectflJJ] to episodic crack propogationOl and the fractal pat- 
terns arising from invasion percolation[3-7]. Drying of col- 
loidal suspensions is also important technologically: paints 
and other coatings depend on colloidal particles for many of 
their key properties, many ceramics go through a stage of par- 
ticle drying, and cosmetics often exploit the unique proper- 
ties of colloidal-scale particles, particularly for such beneficial 
properties as screening the harmful effects of the sun. How- 
ever, for many of these technological applications, the col- 
loidal particles are but one of many different components, and 
drying of the colloids is accompanied by many other phase 
changes. While these mixtures can become highly complex, a 
simpler, yet still rich system that embodies many of the com- 
plex phenomena of these technological suspensions is a mix- 
ture of immiscible fluids with a colloidal suspension; a simple 
example is a mixture of an emulsion and colloidal particles. 
The behavior of the emulsion embodies many of the archety- 
pal phenomena of such systems, while still remaining suffi- 
ciently tractable to enable it to be fully understood. However, 
emulsions themselves typically scatter light significantly, and 
when mixed with a colloidal suspension, this scattering is only 
enhanced. As a result, it is very difficult to image this mixture, 
precluding optical studies of its behavior, and knowledge of 
the actual behavior is woefully missing. 

In this Letter, we explore the drying of mixtures of aqueous 
emulsion droplets and spherical colloidal particles with con- 
focal microscopy, which allows us to resolve the full 3D struc- 
ture of these mixtures and their temporal dynamics. We find 
that the particles first jam into a solidified pack, throughout 
which emulsion drops are dispersed; a front of air then passes 
through the entire system. When this drying front reaches 
large emulsion droplets, the droplets unexpectedly collapse 
and their internal contents are forced into the pore space be- 



tween the surrounding colloids, driven by an imbalance of 
pressures at the droplets' interfaces with air and with the sol- 
vent. By contrast, small droplets are deformed by the dry- 
ing front, yet remain intact without bursting. By coupling the 
Laplace pressure with Darcy's law for flow through a porous 
medium, we predict the duration of large droplet invasion, and 
show that the threshold size between bursting and deformation 
is comparable to the size of the colloidal particles. We use this 
technique to create novel hierarchical materials. 

We suspend sterically-stabilized colloidal spheres of poly- 
methylmethacrylate (PMMA) with radius r p — 1 fjm in dec- 
ahydronaphthalene (DHN). Separately, we create an emulsion 
of an aqueous phase, comprising equal volumes of water and 
glycerol, and PGPR-90 surfactant, required for droplet sta- 
bility, in the nonpolar DHN; our homogenizer creates poly- 
disperse droplets ranging from microns to tens of microns. 
We combine particle suspension and emulsion to create our 
particle-droplet mixtures. These particular components en- 
sure that the refractive indices of the particles, droplets and 
background solvent are all sufficiently matched that we can 
image the entire bulk of the 3D structure with confocal flu- 
orescence microscopy, with single-particle resolution!!]. To 
distinguish particles from droplets, we use different dyes: par- 
ticles are dyed with nitrobenzoxadiazole (NBD) and appear 
green; droplets are dyed with rhodamine-B and appear red. 

We deposit the particle-droplet mixture on a clean glass 
coverslip, and image with an inverted confocal microscope 
(Leica SP5). The droplets are coated immediately by colloids, 
reminiscent of a pickering emulsionfl fioll . As in the case of 
particle-only systems, evaporation proceeds in two stagesfl. 
First, particles and droplets are driven toward the edge of the 
drying sample, where they jam, analogous to the coffee-ring 
effect in particle-only systems[l]. Then, a drying air front 
invades the jammed system and displaces the DHN, which 
ultimately evaporates completely; we do not observe cracks 
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FIG. 1. Structure and drying of droplet-particle mixtures, (a) Recon- 
struction of a typical sample from confocal microscope images, with 
dimensions 91 x 91 x 30 fim 3 . Polydisperse red spheres are emulsion 
droplets; monodisperse green spheres are PMMA particles, shown at 
half size for clarity, (b)-(d) Two-dimensional confocal microscope 
images in x-y plane, within the bulk of the sample, showing the in- 
vasion of the drying front, (b) Onset of air invasion, (c) As the dry- 
ing front moves through, emulsion droplets turn black first, followed 
by (d) air invading the particle regions, (e)-(g) Bright-field images 
showing droplet invasion, where air appears black and the solvent is 
transparent. Air (e) approaches, (f) contacts, and (g) rapidly evacu- 
ates the droplet. 

during drying, as expected for particles of this size^. We use 
confocal microscopy to observe the particle-droplet mixture 
after jamming, as the drying front passes through, which al- 
lows us to determine the exact size and position of each emul- 
sion droplet and colloidal particle [81]. A 3D reconstruction of 
a typical jammed mixture before drying is shown in Fig. 1(a), 
with droplets and particles shown in red and green, respec- 
tively. Collecting this full 3D data takes several seconds, far 
too slow to observe the rapid dynamics that occur during air 
invasion; instead, to capture these dynamics, we fix the focal 
plane in the bulk of the sample and collect high-speed 2D im- 
ages every 70 ms. By imaging the same regions in 3D first, 
then rapidly in 2D as the air invades, we can observe the dry- 
ing process with good spatial and temporal resolution. 

In the jammed configurations before air invasion, the emul- 
sion droplets are distributed throughout the colloidal particles, 
as shown in Fig. 1(a). As the sample dries, more and more par- 
ticle regions turn black; the invading air displaces the solvent 
around the particles, destroying the refractive index match, as 
illustrated in Figs. l(b)-(d). The droplets also turn black as 
drying proceeds, but whether they are replaced by air or sol- 




FIG. 2. (a)-(c) Confocal microscope images of the air invasion of a 
large droplet viewed at a fixed x-y plane, where air forces the con- 
tents of the droplet into the surrounding colloids. Before air invasion, 
tiny projections of the droplet extend into the surrounding pores, vis- 
ible in (a), (d) Schematic in the x-z plane showing the pressure dis- 
tribution before air (white) touches the droplet. Enlargement of the 
(e) air-solvent and (f) droplet-solvent interfaces, (g) Schematic of the 
pressure distribution during evacuation. Inside the droplet, the region 
next to air-droplet interface has a high pressure P2; the region next 
to solvent-droplet interface has a low pressure Pi. The difference 
between these pressures drives the flow of droplet contents into the 
surrounding particles. 

vent cannot be resolved with fluorescence, where both appear 
black. We therefore observe droplet invasion with bright-field 
transmission microscopy, where the solvent remains transpar- 
ent but air appears black. We observe that droplets are re- 
placed by material that appears black in bright-field, which 
therefore must be air, as shown in Figs. l(e)-(g). For a more 
detailed understanding of droplet behavior, we observe large 
isolated emulsion droplets, with radii i?d greater than a few 
microns, during drying. Because it is energetically unfavor- 
able to displace the pre-existing organic solvent, which wets 
the PMMA particles, with the aqueous droplet fluid that does 
not, we should expect that the aqueous fluid remains in the 
droplet and slowly evaporates upon contact with air, with very 
little fluid motion. Contrary to this expectation, however, we 
instead observe significant fluid motion: air rapidly forces the 
contents of the droplet into the pore space between the sur- 
rounding particles, leaving an empty, spherical void, as shown 
in Figs. 2(a)-(c). This flow can not occur spontaneously; in- 
stead, a strong driving force must exist. 

To understand this driving force, we analyze the distribution 



of pressures inside and around an isolated droplet. The total 
pressure near an interface is determined by a combination of 
the external pressure and the interface's Laplace pressure. The 
magnitude of the Laplace pressure is usually estimated as the 
ratio of the interfacial tension to interface's curvature radius; 
its sign positive for convex interfaces, negative for concave. 
Before air invasion, the droplet is surrounded by particles and 
solvent, as shown in Fig. 2(d). The air protrudes into the sol- 
vent, making a concave profile with a negative Laplace pres- 
sure, as shown in Fig. 2(e)|2,]. We estimate P so i v , the pressure 
of the solvent, as the difference between atmospheric pressure 
P a tm, and the Laplace pressure of the air-solvent interface, as 
shown in Fig. 2(e): P so i v = Patra-^airlsoiv/a. where <7 alr | solv 
is the interfacial tension of the air-solvent interface, and a is 
the typical size of the pores between colloidal particles. The 
tiny pore size produces large negative Laplace pressures, and 
hence a low P so iv As a result, the droplet is surrounded by a 
low-pressure environment. 

At the interface between droplet and solvent, the aqueous 
droplet protrudes into the non-polar solvent around the col- 
loids, as shown in Fig. 2(f); the Laplace pressure here is there- 
fore positive. We estimate the pressure inside the droplet as 
Pi = Psoiv + 0drop|soiv/ a > where er dr op| S oiv is the interfacial 
tension of the droplet-solvent interface. 

As soon as air touches the droplet, however, the pressure 
distribution changes dramatically. In particular, part of the 
droplet is now in contact with the much-higher atmospheric 
air pressure, as illustrated in Fig. 2(g). Near the air-droplet 
interface, the pressure inside the drop is estimated as P2 = 
Patm — c air | dl . p/Pid, where we have estimated the radius of 
air-droplet interface by the droplet size Rd- The pressure 
near the droplet-solvent interface remains Pi, as illustrated 
in Fig. 2(g). If P2 > Pi, then the pressure difference will 
force droplet fluid into the surrounding pore space between 
colloidal particles. We estimate this pressure difference: 



AP = P 2 - Pi = 
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where any change in solvent pressure across the droplet is 
negligible. Interestingly, the AP depends not on atmospheric 
pressure, but rather on the competition between the Laplace 
pressures at the various interfaces. We measure the cor- 
responding surface tensions with the pendant drop method: 



r I solv 



26± 2 mN/m, cr a i r |dr 



51 ± 3 mN/m, and 



Cdrop|soiv = 3.8 ± 0.3 mN/m. Because of the presence of 
surfactant at the interface, the droplet-solvent surface tension 
^droplsoiv is so much lower than the other two that its con- 
tribution to the overall pressure difference is negligible. The 
radii that we measure for the large droplets, Rd = 10 — 50 
fim, is orders of magnitude larger than the size of the inter- 
particle pore space, a = 0.36r p = 0.36 /im for random close 
packed particles 11 ill . Consequently, the contribution from a 
pressure drop across the air-droplet interface is also small, and 
the flow is essentially driven by the low pressure in the solvent 
evacuating the large drops, as shown in Figs. 2(a)-(c). Since 
the solvent strongly wets the particles, the menisci of the air 
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FIG. 3. (a) Measurement of droplet evacuation duration time, t cv , 
as a function of droplet radius, Rd, for a typical drying experiment. 
For the measured half decade, the data are consistent with a power- 
law of slope 2, confirming the prediction of eqn.|2] (b) Schematic of 
2D images representing slices through a droplet at different heights, 
(c) Confocal microscope image at the top of the droplet, showing 
the hexagonal crystalline ordering of the particles that suggests (j> ~ 
0.74. (d) Confocal microscope image through the droplet center, 
showing the dense packing around the droplet, consistent with high 
4>. In these images, the brightness of the particles in contact with the 
droplet has been enhanced for clarity. 

evaporating the solvent from the pores creates a low pressure 
in the solvent. It is this, reflected in the first term in Eqn. Q] 
which establishes the large pressure difference that drives the 
flow. We estimate AP = 0.6 atm for these large drops. 

The high time resolution afforded by the rapid collection 
of 2D image sequences in the confocal microscope allows us 
to measure t cv , the evacuation time for air to invade the large 
droplets and force their contents into the surrounding parti- 
cles; i ev varies for droplets of different sizes. We measure the 
size of each droplet from the initial 3D confocal data, then 
quantify t ev using fast 2D images of droplet evacuation. Col- 
lected from deep in the bulk of the sample, our images are 
large enough to contain a number of droplets to obtain good 
statistics, yet are small enough relative to the sample size to 
achieve a uniform sampling environment. We find that the vari- 
ation of t cv with Pd is consistent with a power-law, albeit over 
less than a decade. The exponent is approximately 2, as shown 
by the comparison of the solid line with the data points on the 
log-log plot in Fig. 3(a). 

We estimate this relationship theoretically with a simple 
model based on Darcy's Law for fluid flow in a porous 
medium lll2ll to determine a characteristic evacuation veloc- 
ity, v cv = kVP/ /.i, where k is the permeability of the porous 
medium, and fi is the dynamic viscosity of the fluid. We esti- 
mate VP, the pressure gradient that drives the aqueous fluid 
into the porous medium, as the characteristic pressure differ- 
ence divided by the droplet diameter, VP = AP/2Pd- We 
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FIG. 4. (a)-(c) Behavior of small droplets during drying. Small 
droplets are (a) deformed but not invaded, and eventually (b) coalesce 
to form a continuous phase, (c) The threshold droplet radius, Rth, 
vs. the particle radius, r p . A linear fit through the data points, R t h = 
(0.88 ± 0.16)r p + (0.1 ± 0.1), consistent with our prediction Rth = 
0.83r p . (d)-(f) Hierarchical porous material created from drying. 
Structure made from Rd = 20 /im droplets and r p = 1 /im particles 
(d) before air invasion, and (e) after invasion, (f) Enlargement of the 
compact region marked by the box in (e). The particles appear non- 
spherical because of the difference in their refractive index relative 
to air. The system contains pores of two characteristic sizes: voids 
from droplets (20 /im) and inter-particle pores (0.5 /im). 



then estimate t cv = 2Ra/v cv , yielding: 
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The model predicts that t cv data follow an dependence; 
indeed, the experimental data closely conform to this par- 
ticular power-law scaling, as shown on the log-log plot in 
Fig. 3(a). Moreover, we can further test the model by es- 
timating the prefactor A/j,/(kAP) in Eqn. [2] We estimate 
K using Kozeny-Carman equation for flow through a porous 



medium lll3ll . n = r 2 (l — </>) 3 / (45</> 2 ), where the pore size is 
determined by the particle volume fraction qb. From the 3D 
particle positions in the bulk, we measure <f> = 0.63 ± 0.03, 
consistent with random close-packing; however, around the 
large droplets, the particles are hexagonally close-packed lfiol . 
and the next layer tightly packs the interstices, as illustrated in 
Figs. 3(b)-(d). We think these well-ordered particles predom- 
inantly determine the pore size through which the droplets 
evacuate. We therefore estimate the permeability using the 
higher <fi — 0.74. We measure fj, — 8.9 ± 0.1 mPa-s for 
the H20/glyercol mixture, yielding our rough estimate of 
4/i/(kAP) ps 0.8 x 10 9 s/m 2 ; this value is of the same or- 
der as the value of t ev /R\ = 1.2 x 10 9 s/m 2 from the fit to 
the experimental data in Fig. 3(a), and provides support that 
our model correctly captures the proper physics. 

Our model explains the air invasion of large droplets, whose 
contents are forced into the surrounding pore space. It also 
implies a completely different behavior for sufficiently small 
droplets: as R& decreases, AP must also decrease, as can be 



seen from Eqn. (fTJ. Eventually, for sufficiently small drop 
size, AP will become zero and will no longer drive any flow; 
the high Laplace pressure of these small droplets makes them 
so stiff that they can not be invaded by air. Indeed, when we 
observe small droplets prepared by ultrasonic homogeniza- 
tion, we find that the small droplets are not invaded by air 
but instead are simply deformed and pushed into pore space 
between the surrounding particles, as shown in Figs. 4(a)-(b). 

We estimate the threshold radius P t h, beneath which 
droplets will not be invaded by air by using the previous val- 
ues for the parameters in eqn. (Q~|i and solving for AP = 0; 
our model predicts that P t h = 2.30a = 0.83r p . To test this 
prediction, we measure P t h for droplets mixed with parti- 
cles of several different radii, varying r p by more than half 
an order of magnitude. We find the dependence R t h — 
(0.88 ± 0.16) * r p + (0.1 ± 0.1), in excellent agreement with 
the prediction of our model, as shown in Fig. 4c. 

We observe two qualitatively different behaviors: droplets 
that evacuate and collapse, creating large voids; and droplets 
that remain intact during drying, yielding void-free particle 
packs. We use thisdichotomy to produce hierarchical materi- 
als with several different controllable length scales, by varying 
droplet and particle sizes. One such structure, using droplets 
with i?d = 20 /im and particles with r p = 1 /im, is shown in 
Figs. 4(d)-(f). The resulting hierarchical porous material has 
voids of two length scales: 20 /im, from droplets, and 0.5 /an, 
from the pore space between particles. Hierarchical materials 
may be useful in making low-density porous materials or to 
mimic hierarchical natural structures! 15, our evacuation 
results demonstrate drying as a general low-energy method to 
drive desired materials into a porous medium. 
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